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ABSTRACT 

We present a catalog of 792 DR5 SDSS quasars with optical spectra that 
have been observed serendipitously in the X-rays with the XMM-Newton. These 
quasars cover a redshift range of z = 0.11 - 5.41 and a magnitude range of i 
= 15.3-20.7. Substantial numbers of radio-loud (70) and broad absorption line 
(51) quasars exist within this sample. Significant X-ray detections at > 2a 
account for 87% of the sample (685 quasars), and 473 quasars are detected at 
> 6a, sufficient to allow X-ray spectral fits. For detected sources, ~ 60% have 
X-ray fiuxes between F2_iofceV = 1 - 10 x lO^'^^ ergs cm~^ s~^. We fit a single 
power-law, a fixed power-law with intrinsic absorption left free to vary, and an 
absorbed power-law model to all quasars with X-ray S/N > 6, resulting in a 
weighted mean photon index F = 1.91±0.08, with an intrinsic dispersion a^ 
= 0.38. For the 55 sources (11.6%) that prefer intrinsic absorption, we find a 
weighted mean = 1.5±0.3 x 10^^ cm~^. We find that F correlates significantly 
with optical color, A{g — i), the optical-to- X-ray spectral index (aox) and the X- 
ray luminosity. While the first two correlations can be explained as artefacts of 
undetected intrinsic absorption, the correlation between F and X-ray luminosity 
appears to be a real physical correlation, indicating a pivot in the X-ray slope. 

Subject headings: Galaxies: AGN — accretion disks — X-rays: general 

1. Introduction 

Catalogs are indispensable in performing statistical studies of quasar properties. The 
known correlations between optical and X-ray properties, discussed further below, imply a 
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conn ection between the accretion disk posited to feed the central black hole (IShakura &: Sunyaev 
19731 ) and the ho t Compton-scatterinK cor ona posited to lie in some unknown geometry 
around the disk (IHaardt fc Maraschil Il99ll ) . While quasars were first discovered in early 
radio surveys (e.g. 3C, SCR, PKS, 4C, AO), most quasar surveys since then have been con- 
ducted in the optical/UV. Early quasar selection techniques included UV-excess selection 
and emission line searches using slitless spectroscopy. 

UV-excess selection utilizes the Big Blue Bump that dominates the optical/UV spectrum 
to distinguish quasars from stars. The Big Blue Bump is normally att ributed to a multi- 



temperature accreti on disk ([Shields et a 



19781 : iMalkan fc Sargentlll982f ). UV-excess based 



surveys include the iBraccesi et al.l (Il970l ) catalog, which contains 175 objects with U - B < 
-0.42, and the Palomar-Green (PG) Bright Quasar Survey (BQS. ISchmidt fc Greenlll983l ). 
a sample of 92 objects selected with U - B < -0.44. Slitless spectroscopy, with prisms or 
grisms, obtains a large number of low-resolution spectra for a single field. For example, 
the objective prism on the Curtis Schmidt tel escope at CTIO found 174 corifirmed quasars, 
ranging in redshift from z = 0.1toz = 3.3 (jOsmer fc SmithI Il980l : lOsmerl Il98ll ). Slitless 
spectroscopy often incorporated U V-excess selection as well; the primary example of this 
technique is the Markarian survey (jMarkariaru 119671 : iMarkarian et al.lll98ll ). which searched 
for galaxies with unusually blue continua using a grism. The Large Bright Quasar Survey 
(LBQS; Hewett et al. 1995) used a combination of color selection and the presence of emission 
features on objective prism plates to obtain a homogenous sample of 1055 quasars spanning 
a wide redshift range (0.2 < z < 3.4). 

However, both of these techniques suffer from serious biases. While slitless spectroscopy 
has a high selection efficiency, even at higher redshifts, it is biased against quasars with 
weak emission lines, and cannot reach as faint a flux limit. The UV-excess selection method 
is biased against "red" sources, where red colors may be due to high- redshift (the Lya 
line enters the spectrum at z ~ 2), dust-reddening, signiflcant host galaxy contribution, 
or intrinsically red emission mechanisms. With the advent of the 'UK Schmidt Survey' 
(jWarren et al.lll99ll ). the Two-Degree Field (2DF: ICroom et al.ll200ll ) and the Sloan Digital 
Sky Survey (SDSS; York et al. 2000) quasar catalogs, multicolor selection techniques that 
used up to 5 photometric bands were introduced that could select red quasars in addition 
to blue ones, provided the sources were within the survey flux limit. The 5th Data Release 
(DR5) SDSS quasar catalog has surpassed all previous o ptical surveys by prov iding high 
quality photometry and spectroscopy for 77,429 quasars (ISchneider et al.l 120071 ) spanning 
redshifts from z = 0.08 to z = 5.41. 



As quasars emit strongly over eight decades of the spectrum (e.g. lElvis et al.l Il994j ). 
multiwavelength surveys are necessary to relate the optical/UV accretion emission to other 
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components, notably the non-thermal emission seen in the X-rays. However, X-ray spectra 
are time-consuming and expensive to obtain for large samples. For this reason, previous 
studies of optical and X-ray correlations consist largely of two types: (1) Small samples (N ~ 
20-50) with detailed X-ray spectral analys is have been compiled by observiiig sub-samples of 



optical surveys with X-ray telescop es fe.g . iLaor et al. 



1997 



20051 : IShemmer et all l200d 120081 1. The lAkvlas et all (12004 ) XMM-Newton/2(lY survey is 



Elvis et al. 1994: PiconceUi et al. 



larger, with 96 2QZ quasars observed in wide field (2.5 deg^), shallow (2-10 ks per pointing, 
/ (0.5-8 keV) ~ 10"^^ erg cm~2 s^^) XMM observations. (2) Still larger samples (N ~ 200- 



studies have to ass ume an X-ray spectral slope (e.g. 
Steffen et allbooeh . 



Viffnali et al. 


2003; 


Strateva et al. 


2005 



The ever-expanding archive of X-ray observations now provides a less costly and time- 
consuming rnethod of obtaining X-ray spectra for large numbers of optically-selected quasars. 
Kelly et al.l (120071) cross-corre lated the SDSS DR3 quasar catalog with archival Chandra 
observations (IWeisskopil 1 19911 ) to obtain 174 quasars, 44 of which have sufficient counts to 
fit an absorbed power-l aw with F and as free parameters. This sample was extended to 
a total of 318 qua sars (IKellv et al.l 120081) by adding 149 RQ objects from an SDSS-ROSAT 
cross-correlation ( Strateva et al. 20051 ). As only single power-law models were fit to the 



ROSAT spectra with sufficient counts, 153 of 318 sources have X-ray spectra slopes, but 
no fits for intrinsic absorption were made. Archival Chandra observations have also been 
used to identify an X-ray-selected AGN population in the Chandra Multiwavelength Project 
( 1 Green et al. I l2003h . In addition, a large sample (N = 1135) of optically sel ected SDSS quasar s 
with photometric redshifts are detected in the X-rays in Chandra fields (I Green et al.l |2 0081 ). 
of which 156 have sufficient counts (> 200) to fit an absorbed power-law with F and N// as 
free parameters. 

XMM-Newton is a good choice for cross-correlating with optical catalogs due to its 
large field of view and large effective area. Images a re made with the three European 
Photon Imaging Cameras (EPIC): MOS-1 and M OS-2 (iTurner et al.ll2001l) . each of which 
has a 33' x 33' field of view, and the pn camera (IStrueder et al.ll200ll ). which has a 27.5' 
X 27.5' field of viewfl. XMM's large effective area (922 cm^ for MOS and 1,227 cm^ for 
PN at 1 keV)^ results in higher signal-to-noise X-ray spectra than with Chandra for bright, 
non-background-limited source^. The EPIC CCDs have good spectral resolution (E/AE ~ 
20 - 50 for both MOS and PN) over the 0.5 - 10 keV band. 



^http: / /imagine. gsfc.nasa.gov/docs/sats_n_data/missions/xmm. html 
^http://heasarc. nasa.gov/docs/xmm/uhb/node86.html 
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Archival XMM-Newton observations overlapped with ~1% of the SDSS DR5 coverage 
as of Feb. 2 007. The SP S S and XMM-Newton archives are well -matched in sensitivity, as 



discu ssed in lYoung et al.l (120081 ). The 5th Data Release (DR5, lAdelman-McCarthy et al. 



20071 ) of the SDSS covers a spectroscopic area of 5740 deg^, and contains 90,611 quasars. 
The SDSS photometry in the u, g, r, i and z bands covers 3,250 - 10,000 A, while the 
spectroscopy covers 3,800 - 9,200 A with a spectral resolving power of ~ 2,000. The SDSS 
is ~95% co mplete for point sources to a limiting magnitu de i = 19.1, corrected for Galactic 
reddening (IRichards et al.l 120021 : 1 Vanden Berk et al.ll2005l ). 



A preliminary cross-correlation of DRl SDSS quasars with t he XMM-Newton p ublic 
archive yielded 55 objects with expo sure times greater th an 20 ks (IRisaliti et al.ll2005l ). Of 
these, 35 yielded good X-ray spectra. IRisaliti et al.l (120051 ) estimated that a cross- correlation 
of a final SDSS data release with the ever-growing XMM-Newton archive would yield ~ 1000 
quasars, of which ~80% would have good X-ray spectra. 

In this paper, we cross-correlate the SDSS DR5 quasar catalog with archival XMM- 
Newton observations to obtain a large (N~800) sample of quasars with X-ray detections, 
almost 500 of which have good optical and X-ray spectral data. Below, we outline two 
immediate goals for the SDSS/XMM-Newton Quasar Survey: (1) to conduct large, statis- 
tical studies to understand the physical basis behind optical/X-ray correlations and (2) to 
investigate interesting sub-populations of quasars. 



1.1. 



Optical/X-ray Correlations 



The relations between optical and X-ray continuum and spectral properties promise to 
reveal clues about the disk-corona structure of quasars. The large sample provided by the 
SDSS/XMM-Newton Quasar Survey allow two correlations to be investigated. The first of 
these is the controversial a„r— ^ ° relation: many studies have found that the spectral index 

2500A ^ 

from 2500 A to 2 keV, defined as a„o 



\og(L2kev/ L2500)/log(z/2fcfi 



with the log of the monochromatic luminosity at 2500 A, l^^^^^ (Tananbaum et al 

Kriss Sz CanlzaresI 



Zamorani et al. 



1981 



Anderson &: Margonll 19871; 



Avni Sz Tananbaum 



1982 



/ ^25nn) , anti-correlates 



1985 



19941: IPickering et al.lll993; 



1979 



Tananbaum et al. 



1986 



1995 



Vignali et al. 



Wilkes et al 



2003l: lstrateva et al.lbood: 



20071 ). However, some studies (IBechtold et al. 



Shen et al 



2003 



Avni. Worrall &: Margon 



20061: Isteffen et al.ll2006l : ljust et al. 



tion to be between a or and redshift, while other studies (Yuan et al.l 



Kelly et al.ll2007) find the primary rela- 



19981 : iTang et al.ll2007h 



find that the correlation may be induced by selection effects. 



No physical basis for the aox — ^gsooi ^^^^^^^^ been proposed, and the relation 
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itself provides little guidance. In part, this is because previous studies have largely used the 
traditional, observationally convenient, but physically arbitrary, endpoints of 2500 A and 2 
keV. They also assume an X-ray photon index (F ~ 2, where T = -a + 1 for Fj, oc z/°) to 
obtain the X-ray flux at 2 keV. A systematic study with measured optical and X-ray spectra 
would enable an investigation of the relation at different frequencies than those traditionally 
used, hopefully revealing clues about the relation's physical underpinnings (Young et al. 
2009b, in prep.). 

The second c orrelation is the posi t ive re lation between F and the normalized accre- 
tion rate, L/LEdd Jshemmer et al.lbood . boosh . While early studies focu sed on the relation 



between F and full- width half-max imum of the H/3 line, FWHM(H/3) (IBoUer et al.l Il996 



Laor et al.lll997l : iBrandt et al.lll997l ). Shemmer et al. (2006, 2008) have broken the degener- 
acy between FWHM(H/?) and L/LEdd by including highly luminous sources in order to show 
that F depends primarily on L/LEdd- There have been suggestions that this dependence may 
be due to the disk emitting more and softer p hotons as accretion rates increase, leading to 
more efficient Compton cooling in the corona (jLaoiil2000l : iKawaguchi et al.ll200ll ). 



Until now, studies of the T-L/LEdd relation have consisted of small samples (N ~ 40). 
The SDSS/XMM-Newton quasar survey can increase this sample size by an order of magni- 
tude, leading to a better defined relation (Risaliti et al. 2009, in prep.). 



1.2. Quasar Sub- Populations 



Interesting quasar sub-populations can be readily investigated due to the large number 
of sources in the SD SS/XMM-Newton Qu asar Survey. For example, red quasars make up 6% 
of the SDSS sa mple (IRichards et al.ll2003f). Their steep op tical slopes have been attributed to 



dust-reddening (IRichards et al. 



2003 



Hopkins et al.ll2004l ). th ough observations o f individual 



objects s uggest that some slo pes may be intrinsically steep ( Risaliti et al. 20031 : Hall et al. 



20061 ). In lYoung et al.l (120081 ) . we studied 17 quasars with extreme red colors {g — r > 0.5) 
and moderate redshifts (1 < z < 2). By using X-ray observations in conjunction with optical 
spectra, we constrained the amount of intrinsic absorption in each source, thereby allowing 
the separation of intrinsically red from dust-reddened optical continua. We find that almost 
half (7 of 17) of the quasars can be classified as probable 'intrinsically red' objects. These 
quasars have unusually broad Mgll emission lines (<FWHM>=10,500 km s~^), fiat but 
unabsorbed X-ray spectra (< F > = 1.66±0.08), and low accretion rates {M / MEdd ~ 0.01). 

Other interesting sub-populations for future investigations include broad absorption line 
(BAL), Type 2, and radio-loud (RL) quasars. 
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In this paper, we describe the SDSS quasar selection and the method with which we 
match sources to XMM-Newton observations in §2. X-ray data reduction is described in 
§3 and the resulting sample and correlations are discussed in §4. We assume a standard 
cosmology througho ut the paper, where Hq = 70 km s~^ Mpc~\ Qm = 0.3, and = 0.7 
JSpergel et al.ll2003h . 



2. Data 



2.1. SDSS Quasar Selection 



As the DR5 quasar catalog (ISchneider et al.l 120071 ) was not yet available at the time of 
selection, we selected quasars with optical spectra directly from the DR5 SDSS databas^ by 
choosing SpecClass = 3 (QSO) or 4 (QSO with z > 2.3, whose redshift has been confirmed 
using a Lya estimator). These quasars were selected for spectroscopic follow-up by the SDSS 
primarily due to their photometric col ors, although some quasars were selected because they 
have a match in the FIRST survey (jWhite et al.l 119971 ). No X-ray selection is involved. 
Target selection efficiency, i.e. the percentage of sources spectroscopically ob served that are 
confirmed as quasars, is 66% jRichards etaPbooi IVanden Berk et al.l 120051 ). The selection 
method is described in detail by 



Richards et al 



and is summarized briefly below. 



The vast majority (~95%) of SDSS quasar candidates are detected using a multicolor 
selection technique. Quasar candidates are defined to be any object at least 4cr away from 
the stellar locus, which is defined in the [u — g,g — r,r — i,i — z) color space. In addition, 
special color-color regions are defined to specifically include or exclude quasar candidates. 
Inclusion regions include quasar candidates from 2.5 < z < 3, even if their colors cross the 
stellar locus. Exclusion regions prevent contamination due to white dwarfs, A stars and M 
star-white dwarf pairs. 

For both radio and color-selected quasar candidates, magnitude limits are applied. 
Quasar candidates brighter than i = 15 are rejected for spectroscopic follow-up because 
bright sources can contaminate the spectra of objects in adjacent fibers in the SDSS spec- 
tograph. Radio and low-redshift color-selected candidates fainter than i = 19.1, which have 
a high number density on the sky, are also rejected due to the limited number of optical 
fibers for follow-up spectroscopy. Since high redshift (z > 3) quasar candidates have a lower 
surface density on the sky, a fainter cut-off magnitude i = 20.2 is applied to these objects. 



■^http:/ / www.sdss.org/dr5/access/index.html 



- 7- 



SDSS selected a small number of quasar candidates (~5%) by matching point sources 
with the position of a radio detection fr om the FIRST survey to within 2". While the 
DR5 quasar catalog (ISchneider et al.ll2007l ) includes sour ces selected by rnatchi ng to ROSAT 
detections, the initial SDSS quasar selection outlined in [Richards et al.l (120021 ) does not use 
X-ray detection criterion. 

While SDSS radio selection requires that a source be point-like, their color selection 
includes extended sources as well, in order to include low-redshift AGN, such as Seyfert 
galaxies. However, in addition to having colors distinct from the stellar locus, extended 
sources must also have colors distinct from the main galaxy distribution. (The main galaxy 
distribution overlaps the stellar locus; however, a galaxy can be a clear outlier from the stellar 
locus both due to the shape of the stellar locus and because the stellar locus is determined 
primarily from F and M stars that dominate the Galactic stellar density at high latitudes.) 
Simple color cuts are applied to distinguish extended sources from galaxies, rather than an 
additional multicolor selection. 



The DR5 quasar catalog (ISchneider et al.ll2007l ) additionally requires: luminosities brighter 



than Mj = -22.0; at least one emission line with a FWHM greater than 1000 km s~^ or in- 
teresting/complex absorption features; magnitudes fainter than i = 15.0; and highly reliable 
redshifts. This results in a catalog of 77,429 quasars. We used only those quasars that are 
also in the DR5 Quasar Catalog for further analysis. Of the 92 sources we reject from our 
sample, ~40% are Type 2 quasars. These will be the subject of a later study. 

We matched The selected SDSS quasars with the FIRST survey Jwhite et al.l [l997llR 
using a search radius of 3" in order to calculate their radio-loudness (Rl = ^^ghz / F4400, 
Kellerman et al. 1989)). A source is taken to be radio loud (RL) if R^ > 10. A power-law is 
interpolated between the optical magnitudes to get Fa(4400A), and the 1.4 GHz radio flux 
is obtained from FIRST survey detections, which are extrapolated to 5 GHz using a radio 
power-law ur = -0.8. All the quasars lie in the area covered by the FIRST survey, so if there 
is no detection, we use the 5a upper-limit on the 1.4 GHz radio flux to extrapolate to 5 GHz. 
For 40 quasars, the upper-limit is too high to determine whether the source is radio-loud. Of 
the remaining quasars, 70 (9.3%) are radio-loud and 682 (90.7%) are radio-quiet, a typical 
ratio (e.g. Kembhavi & Narlikar 1999, pp. 256-263). 

The SDSS c olor-color selection is effective at finding a variety of Br oad Absorption Line 
(BAL) quasars JSchneider et al.ll2007h . The official DR5 BAL catalog Joibson et allboogh 
had not yet been published during the writin g of this paper, so we used an incomplete list 
of ~ 4200 BALs in the DR5 Quasar Catalog (jShen et al.ll2008al ). selected using traditional 



"'http:/ /sundog. stsci.edu/ 
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criteria (IWeymann et al.lll99ll ). Using IShen et al.l (l2008al ). we find 52 BALs in the SDSS- 
XMM Quasar Survey. Of tfiese, 15 BALs liave fiigli enougli X-ray signal-to-noise to obtain 
X-ray spectra (see §3). Since BALs are normally identified by CIV absorption features, 
which are visible only above a redshift z = 1.6 in SDSS spectra, BALs are unlikely to be 
identified in -56% of the SDSS-XMM quasars. BALs make up 14. 5% of the SDSS/X MM 
quasars with z > 1.6, which is in line with the parent population in lShen et al.l (j2008al ). 



In this paper, we measure the optical color of a quasar by its relative color (IRichards et al. 



20031 ). Relative colors compare a quasar's measured colors with the median colors in its red- 
shift bin, where redshift bin sizes are 0.1 in redshift, so that ^{g — i) = {g — i) - <{g — i)>z- 
The use of relative colors corrects for the effect of typical emission lines on the photometry 
in a particular band. The relative [g — i) colors of the SDSS quasars match a Gaussian 
distribution on the bl ue side but require the addition of a tail on the red side [see Fig. 3 in 
Richards et aP J2003h ]. 



2.2. Matching with the XMM-Newton Archive 

We matched the SDSS quasars with the XMM-Newton archive from February 2007, 
choosing only those quasars that fell within 14' of XMM observation field centers (typically 
2-3 quasars per field). As part of the extraction process, described below, a source region is 
defined around each set of SDSS quasar coordinates. Depending on the S/N of the source, 
the extraction radius can range from 10 - 85", with a typical radius of 19". Low S/N objects 
are extracted with smaller radii to minimize the effect of high background levels, while larger 
radii for high S/N objects allow for an increased encircled energy fraction in the presence 
of relatively low background levels. The large extraction radii take into account the XMM 
point spread function (PSF), which is characterized by the radius at which 90% of the total 
energy is encircled. This radius increases from 48" (MOS) and 51.5" (PN) at 0' off-axis 
angle to 52.5" (MOS) and 66" (PN) at 12' off-axis anglfl The extraction radii also take 
positional accuracy into account: SDSS positional errors are negligible (0.1" at the survey 
limit of r=22 for typical seeing. Pier et al. 2003), while the XMM positi onal accuracy is 3 " 



at 3(7 for offset angles < < 10 arcmin and 6" for 5 < < 10 arcmin (jPierre et al.l 120071 ). 



Multiple X-ray observations exist for 265 sources. In these cases, all observations were 
retrieved and reduced, but only the observation with the longest exposure time was used for 
further analysis. To avoid biases, we did not select the highest S/N observations, although 
in ~80% cases, the two selections would be effectively the same. 



^http: / /xmm. vilspa.esa.es/external/xnim_user .support/documentation /uhb/nodel 7.html 
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The chance of including an unrelated random source within the extraction region is small 
but non-zero. Within a ~14' radius field of view, and using the average extraction radius, 

ce there are ~70 sources in a typical 



[ "beams" in j 


m observation. 


(Watson et al. 


2008, which has 



that in this paper), this results in a 3.6% chance of extracting a random source rather than 
the SDSS-selected quasar. Since a random source is likely to be faint in the X-rays, any 
contamination is only significant for sources under 100 net counts. We have 792 unique 
sources, of which 390 have less than 100 net counts, so ~14 (2% of the SDSS/XMM quasar 
sample) have significant contamination from an unrelated 'interloper' source. 



2.3. X-ray Data Reduction 

The 582 X-ray observations were processed using the XMM-Newton Science Analysis 
System, SAS v7.02§. We reprocessed the events to ensure that each observation has the 
same, up-to-date calibration, and then filtered the observations to remove time intervals of 
flaring high-energy background events using the standard cut-off of 0.35 cts s^^ for the MOS 
cameras and 1.0 ct s'^ for the PN cameralll. Source and background regions were defined 
in a semi-automatic process. The SAS task eregionanalyse was used to optimize the source 
extraction radius for signal-to-noise. Most radii include at least 80% of the source counts. 
Background regions were defined by eye, avoiding obvious X-ray sources and chip edges. 
These regions were typicafly a circle of radius 2000-2500 pixels (100-125"), selected to lie at 
the same off-axis angle as the source and as close to the source as possible without overlapping 
the source extraction region. Once source and background regions were defined for every 
SDSS quasar in an observation, spectra were extracted for a total of 1380 non-unique quasars 
in 582 observations. 

To check for biases in the data reduction. Figures la — b show the net counts and the 
X-ray photon index (F) plotted against the extraction radius. Figures Ic — d plot the off-axis 
angle and the X-ray S/N against F. In Fig. la, the net counts are expected to correlate 
with the extraction radius, since the extraction radius will increase to larger encircled energy 
fraction for sources that stand higher above the background. The lack of correlations in Fig. 
Ib-c show that the encircled energy correction takes the extraction radius and off-axis angle 
into account correctly. Fig Id shows that objects with flat F are primarily found among low 
S/N objects, where absorption may be undetected in a spectral flt. 



^http:/ /xmm.esac.esa.int/sas 

^http: / / xmm.vilspa.esa.es / sas/7. 1.0/ documentation/threads / 
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Where possible, observations were processed for all three XMM EPIC CCDs. In ~40% 
of the observations, a source lies in a bad region in one or two of the three cameras, either 
in a strip between two chips or, because the MOS and PN cameras have different shapes, 
outside the field of view in one of the cameras. In these cases, we use the remaining images 
from the other cameras for analysis. 

Table [1] contains observational data for each quasar in the DR5 SDSS/XMM- Newton 
Survey: the SDSS name, XMM observation ID, redshift. Galactic column density in the 
direction of the source. X-ray signal-to-noise, observation exposure time, off-axis angle, net 
source counts, background counts (from a background region that is scaled to the area of 
the source region), and two flags indicating if a source is RL or BAL. 

Figure [2] summarizes the survey characteristics. The X-ray exposure times (Fig. 2a) 
range from 1.6 to 294 kiloseconds, though the majority of observations lie between 20 and 
100 kiloseconds. This range in exposure times results in a wide range in sensitivity. While 
most sources have low signal-to-noise, a significant fraction have S/N > 10, where more 
complex models can be fit (Fig. 26). The detection fraction is > 80% until z > 3.5 (Fig. 
2c), but spectral coverage drops off fairly quickly for sources with z > 2 (Fig. 2d). 



3. X-ray Analysis: Spectral Fits 

We made fits to the extracted spectra using the Sherpa packag^ within CIAo|^. For 
each source, the available MOS-I-PN spectra were fit simultaneously over the 0.5 - 10 keV 
band. The observations were fit according to their S/N, with more complicated models being 
applied as S/N increased. All the models included local absorption fixed to the Galactic 
hydrogen column density (NH,gai) at each source location. Values for l^H,gai were taken from 
the N/f tool available at WebPIMM^^, which is based on the 21 cm HI compilation of Dickey 
& Lockman (1990) and Kalberla et al. (2005). 

For the 319 low S/N sources (S/N < 6), we fix a power-law to the weighted mean 
obtained for the high S/N (S/N > 6) quasars in the sample, F ~ 1.9 (§4.1), and allow only 
the normalization to vary in order to obtain the fiux. For the 101 sources with S/N < 2, we 
obtain a 90% upper-limit to the fiux. We use the Cash (1979) statistic, which gives more 
reliable results for low-count sources, to fit sources with S/N < 6. As the background is 



http:/ / cxc.harvard.edu/ sherpa/threads/index.html 
http://cxc.harvard.edu/ciao/ 

http: / /heasarc. gsfc.nasa.gov/Tools / w3pinims.html 
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not subtracted and is instead fit simultaneous ly with the s ource , we apply a background 



model with three components, as described in iLumb et al.l (120021 ) and in the XMM Users 
Handbool0: a power-law (for the extragalactic X-ray spectrum), a broken power-law (for 
the quiescent soft proton spectrum), and two spectral lines (for cosmic-ray interactions with 
the detector). All parameters except for the normalization of the spectral lines were fixed. 
For 13 undetected sources, and 5 detected sources, this fit results in bad or null flux values. 
In these cases, we flag the source as undetected (flag = -1 in column (9) of Table [2]) and we 
list the values as 9.99 (column (8) of same table). 

For the 473 sources with enough S/N to fit spectral parameters (S/N > 6), we use the 
statistic to fit three models: a single power-law (SPL) with no intrinisc absorption, a 
fixed power-law (FPL) where intrinsic absorption is left free to vary, and an intrinsically 
absorbed power-law (APL). The F-test0, which measures the significance of the change in 
as components are added to a model, is used to determine whether the data prefer the 
APL model. To compare the SPL and FPL models, we simply compare the respective 
values, since both models have the same number of parameters. Therefore, the best-fit model 
is the one preferred by the F-test that also has the lowest value. 

For sources with an unacceptable x^ values for all three models, we by default make the 
SPL model the best-fit, but we also list the APL 90% upper- limit on intrinsic absorption. 
We plot the reduced distribution in Figure [31 A later paper will look in more detail at 
sources with bad fits (x^ > 1-2). 

Results for the best-fit models are listed in Table [21 including observed-frame and rest- 
frame fluxes (or 90% upper-limits), ctox, a flag indicating the best-fit X-ray spectral model, 
the photon index, intrinsic absorption (or 90% upper-limit) and the x^ values and degrees 
of freedom for the best-fit model. The best-fit flag indicates which values are listed for each 
best-fit model. For sources that prefer the APL model (flag = 3), both P and Njy are from 
the APL fit. If the FPL model is preferred (flag = 2), the SPL P and FPL N/^ are listed. 
For sources that prefer the SPL model (flag = 1), P is from the SPL fit and Nj^ is the 90% 
upper-limit from the APL fit. 



^^littp: / /xmm. vilspa.esa.es/external/xmm_user_ support/documentation 
^ ^ http : / / cxc. harvard, edu /ciao/ ahelp / ft est .html 
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Results and Discussion 



The SDSS/XMM-Newton Quasar Survey contains 792 sources, 685 of which are detected 
in the X-rays and 473 of which have X-ray spectra. (All have optical spectra.) The catalog 
covers redshifts z = 0.11 - 5.41 and optical magnitudes range from i = 15.3 to i = 20.7. 
Figures 4 shows the survey sensitivity in the optical (a) and X-ray (6) bands, and the 
observed-frame 2-10 keV flux distribution (c). The FWHM of the distribution spans F2-iofeeV 
= (1 - 10) X 10^^^ ergs cm~^ s"'^ and contains ~ 60% of the detected sources. 

Results of standard aox analysis (i.e. using the conventional 2500 A and 2 keV fiducial 
points) are shown in Figure [5l We test for correlations with ap^, using the K endall correlation 
test available in ASURV, a survival statistics package (ILavalley et al.lll992l ). Figure 5a shows 
log L2keV vs. log L2500 with a dotted line of slope unity for reference. The best-fit line, 
obtained from the EM (estimate and maximize) algorithm within ASURV, is flatter than 
unity with a slope of 0.64 ± 0.03. This deviation from unity is clear in Figure 56, which 
shows the aox ~ ^2500! ^o^^^^^^io^' significant at the llAa level in the SDSS/XMM Quasar 
Survey. The solid line is the best-fit EM regression line for our data: 



Or. 



(3.080 ± 0.376) + (-0.153 ± 0.012)log{L 



The slope is consistent within errors with the best-fit regression from ISteffen et al.l (120061 ) 
(dotted line in Fig. 5b). 

Figure 5c plots aox against redshift. We find no correlation between aox a nd redshift, as 



2003; 


Strateva et al. 


2005; 


Risaliti et al. 


2005 


2003; 


Shen et al. 


20061: 


Kellv et al. 


200 


7). I 



20061) but not all flBechtold et al. 



Kelly et al.l (120071 ) apply a more 



sophisticated statistical analysis by allowing for non-linear fits with multiple variables. As 
a result, they find that aox depends on both Z2500 and redshift so that quasars become more 
X-ray loud at low luminosities and higher redshifts. We have applied only linear regressions 
in this paper, and will apply more complex statistics in a later publication. 

The weighted mean of the X-ray spectral slope for the 473 sources with an X-ray S/N > 
6 is < r > = 1.91±0.08, with a standard deviation of 0.40, which is consistent with previous 
results (e.g.. IShen et al.ll2006l ). The typical la error on F is 0.15, resulting in an intrinsic 
dispersion err = 0.37. Figure [6] shows the aox and F distributions. 

We calculate the percentage of sources with significant intrinsic absorption by testing 
whether sources with S/N > 6 prefer the APL model over the SPL model, with F-test 
probability Pf > 0.95 and acceptable xl, resulting in 34 (7.2%) absorbed sources. However, 
this method is biased against low S/N sources, since the APL model has more parameters 
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than the SPL model. Since the FPL and SPL models have the same number of parameters, 
we count sources as absorbed if they prefer either the FPL or APL models with Pi? > 0.95 
and acceptable xl- This method gives 55 sources (11.6%) that are intrinsically absorbed, 
which is comparabl e to the percentag e found among Type 1 quasars in previous surveys 
(iMateos et aD bood : ICreen et all boosl) The XMM-Newton COSMOS survey obtained a 



Mainieri et al. 



20071 ) using a lower confidence threshold (P^ > 0.9) 



higher percentage (20%, 
when applying the F-test. When we redo our method using the same confidence threshold, 
we find 82 absorbed sources (17.4%). 

Nevertheless, the amount of absorption in this sample is only a lower-limit for two 
reasons. First, undetected absorption may still exist, particularly in low S/N sources. For 
example, in the correlation plots (e.g. Figures 7,9 and 11), three sources can be seen with 
r ~ 0.5. These sources do not formally prefer either absorption model, but their X-ray 
spectra are cut off at soft energies, suggesting absorption as the likely cause of flat spectra. 
Second, we cannot test for absorption in sources with S/N < 6, and it is possible that these 
low S/N sources have a higher percentage of absorption. 

For the 55 sources that prefer absorption , we calculate the weighted mean of the intrinsic 
column density: Nh = 1.5±0.3 x 10^^ cm~^. The absorbed population consists of 6 broad 
absorption line objects (BAL), 7 radio-loud (RL), and 42 radio-quiet (RQ), non-BAL objects. 

Table [3] summarizes the weighted means of Oox, P and for the three sub-sets of 
quasar populations: RQ+non-BAL, RL-|-non-BAL, and BAL. (There are four quasars that 
are both RL and BALs.) Because of the small numbers of absorbed spectra in the sub- 
populations, comparing the respective distributions is not meaningful, but Uox and P, 
Kolmogorov-Smirnov (K-S) tests show that the RQ distribution is significantly different from 
the RL distribution, with probabilities Pks < 0.5% that the two samples are drawn from 
the same parent distribution (Fi g. [SI). RL quasars are known to be brighter in the X-rays for 



a given 2500 A luminosity (e.g. IZamorani et al.lll981l ). However, while RL quasars are also 



known to have flatter X-ray slopes than RQ quasars (IWil 



19921 : iReeves et al.lll997l : iReeves fc Turned l2000l : IPage et al 



ces fc 



2005 



Ivis 1987; iWil 



Piconcelli et al 



iams et al. 



2005|), the 



average P value for RL quasars in the SDSS/XMM sample (< P > = 1.85 ± 0.04) is steeper 
than that found in previous samples (< P > ~ 1.5-1.75). Figure 66 shows that RL quasars 
do follow the same Gaussian distribution as RQ quasars for P < 1.9, but for P > 1.9, the 
RL distribution falls off rapidly. Only 7 of 49 RL quasars (14%) have P > 2, compared to 
48% of RQ quasars. 

The P distribution is not significantly different for BAL vs. non-BAL quasars {Prs 
= 75%) once the X-ray spectra are corrected for intrinsic absorption. However, the aox 
distribution of BAL quasars, which includes 37 sources with S/N < 6 that cannot be corrected 
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for absorption, remains significantly different [Prs = 0.2%) from non-BAL quasars. Previous 



studi es liave noted tlie difference in X-ray briglitness for BAL quasars (e.g. iGreen &: Matfiur 
1996h . 



4.1. Correlations with the X-ray Spectral Slope 

The large number of X-ray spectra in the SDSS/XMM Quasar Survey allows us to test 
for correlations with the X-ray photon index (F). Since we use only sources detected with 
a high enough S/N to fit a power-law, we do not use survival statistics and instead use the 
Kendall's rank correlation coefficient to test for correlations. We use only those sources that 
do not prefer either absorption model, and we also require that the SPL model result in a 
reasonably good fit (x^ < 1.2). In addition, we restrict the selection to those sources with 
S/N > 10 in order to reduce the effect of undetected absorption in lower S/N spectra. 

When correlations are significant, we plot a Weighted Least Squares (WLS) regression 
line to take into account the measurement errors in P, which are much larger than the errors 
in the independent variable. RQ quasars show two correlations of P with luminosity and 
optical color with probabilities less than 0.5% that they are due to chance: (1) P vs. log 
L2keV (Pk = 2.4e-5) and (2) P vs. A(g — i) (Pk = 7.5e-4), where Pk is the Kendall rank 
two-sided significance level for each correlation. RQ quasars also show a marginal correlation 
between P and aox (Pk = 1-6%). The P - L2kev, T - A{g — i) and P - aox relations are not 
significant {Pk > 10%) for RL quasars (green stars in Figures 7, 9 and 11). We now discuss 
each of these correlations in turn. 



4.1.1. X-ray Slope vs. X-ray Luminosity 

An anti-correlation exists between the X-ray slope (F) and the 2 keV luminosity (log 
\j2kev)i hereafter hkeVi such that X-ray slopes harden as X-ray luminosity increases (Fig. 
7a). The WLS regression line (for RQ, non-BAL quasars) is: 

F = (6.39 ± 1.3) + (-0.16 ± 0.05)W (1) 

To investigate the possibility of a pivot point in the X-ray spectrum, we next examined 7 
additional correlations between the X-ray slope (F) and the monochromatic X-ray luminosi- 
ties at 0.7, 1, 1.5, 4, 7, 10 and 20 keV. The monochromatic fiuxes are obtained for sources 
preferring the SPL model by normalizing the fit at each energy in turn. The results of these 
correlations are summarized in Figures 7a — d and in Table HI The strength of the correlation 
increases with energy as the slope steepens, so the strongest, steepest correlation is between 
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r and hokev (Fig. 7d). Correlation strength decreases at lower energies, bottoming out at 
1 keV, where the slope is consistent with zero. At 0.7 keV, the slope flips to a positive 
correlation and the strength of the correlation increases again. The flip in the correlation 
slope indicates a pivot in the X-ray spectrum near 1 keV (Fig. 8a). 

For the correlations above, we have used F to extrapolate the rest-frame monochromatic 
luminosities for sources with redshifts out of range of the observed spectrum. To check 
that this extrapolation does not affect the correlations, we perform the fits again, this time 
excluding any sources where the rest-frame luminosity is not in the observed range. The 
results are shown in Table 5, and plotted in Fig. 86. The errors are larger, due to the smaller 
sample sizes and the narrower range of luminosities observed, but the spectrum still pivots 
between 1 and 1.5 keV. At the lowest energies, the slopes become much steeper, possibly 
due to the influence of a soft excess component on the power-law fit. 

The F determined from the SPL model is used to determine the monochromatic flux, 
so it is important to check that the model assumptions do not induce the observed anti- 
correlation. To do this, we re-fit the sources, this time with a power-law fixed to the sample 
mean (F = 1.91) and intrinsic absorption fixed to zero. We then obtain the monochromatic 
flux at 2 keV and compare the F2kev obtained via a fixed power-law to the F2kev obtained 
via the best-fit power-law. We find that changing the model assumptions changes the log 
flux values by 1.4%, which is not enough to explain the observed correlation at 2 keV. 

Even with the selection restricted to sources with S/N > 10, there is still the possibility 
of undetected absorption. For example. Figures 7 show two sources with F ~ 0.5 and another 
four radio-quiet sources with F ~ 1-1.5. The X-ray spectra of these sources show curvature 
in the soft X-rays that, while not significant enough for the sources to prefer an absorption 
model, nevertheless suggests that intrinsic absorption is the likely cause of flat X-ray slopes. 
Therefore, we test for correlations again, this time using sources with S/N > 20 in order to 
minimize the chance of undetected absorption. The significance of the correlation increases, 
with the probability of a chance correlation falling below Pj^ = le-7. However, with many 
low-luminosity sources eliminated, the correlation is biased to a steeper slope. Having shown 
that the correlations do not rely on undetected absorption, we continue to use the correlations 
for sources with S/N > 10 for further analysis. 



m 



The slope o f the F — hkev anti-correlation is equivalent within errors to the slope found 



Green et al.l (120081 ). They also find a correlation between F and the 2 keV luminosity 
based on 156 RQ and RL quasars fit with a power-law plus intrinsic absorption, plus 979 
quasars with lower S/N that were fit with a single power-law and zero in trinsic abso r ption . 



A similar anti-correlation bet ween F and 2 - 10 ke V luminosity was found in lPage et al.l (120051 ) 



for 16 RL quasars. While the iPage et al.l (120051 ) correlation was not found to be significant 
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for RQ quasars, there were only 7 RQ quasars in their sample. However, a previous study by 
Dai et al.l (120041 ) found the opposite correlation between T and /2-iofcey at 98.6% significance 
for a sample of 10 quasars obs erved with Chandra and XMM. The F and extrapolated hkev 
values for the iDai et al.l (12004 ) quasars are plotted in Fig. 7a as open^lack squares, where 
they are consistent with the trend observed in SDSS-XMM quasarsO The SDSS/XMM- 
Newton Quasar Survey increases the sample size of previous studies by factors of 3 - 30, and 
covers ~ 3 decades of X-ray luminosity. 

To explain the observed correlations between F and X-ray luminosity, we must answer 
two questions: 

1. Why does the X-ray slope change with luminosity? 

2. Why does the slope change such that the pivot point is near 1 keV? 

We address two possible answers to the first question. First, an increased hard compo- 
nent at higher X-ray luminosities may explain the trends observed h ere. The hard compo- 



nent may be due to nonthermal emission associated wit h a jet (e.g. 



Wilkes fc Elvid 119871 ) or due to a reflection component (IKrolik Ill999l ). RL quasars do not 



Zamorani et al. 1981 



show a significant correlation between F and Lx except for marginal correlations at high en- 
ergies (7 and 10 keV), but this may be because a jet already dominates the X-ray emission. 
The spectrum of the hard component may flatten the X-ray spectrum by covering up the 
steeper power-law due to inverse Compton scattering, while simultaneously increasing the 
2-10 keV luminosity. 

Alternatively, quasars with low X-ray luminosities may have steeper slopes due to a 
component linked to high accretion rates. As discussed in §1.1, previous studies have found 
a strong correlation between F and the Eddington ratio (Lhoi/^Edd), where steeper sources 
are associated with higher accretion rates. Therefore, sources with high accretion rates (and 
steep X-ray spectra) must be associated with low X-ray luminosity in order to produce the 
T-Lx correlations. However, it is not clear if high accretion rates and steep X-ray spectra 
tend to coincide with low X-ray luminosity. Studies of X-ray binary (XRB) accretion states 
have shown that high accretion states tend to be a ssociated with steep X-ray spec tra and high 
X-ray flux compared to the low accretion state ( Remillard fc McClintocjbood ). However, 



the distinction is not clear in every source - even "high" accretion states can be associated 
with both low and high X-ray flux. Narrow-line Seyfert 1 (NLSl) galaxies are an extreme 



^•^More recently, ISaez et al.l (j2008l ) find the same correlation to be significan t to > 99.5% for 173 bright, 
radio-quiet quasars in the Chandra Deep Fields. However, the trends found in Saez et al. ( 2008h are domi- 
nated by Type 2 quasars, particularly at low X-ray luminosities. As we have excluded Type 2 quasars, the 
two samples are not in conflict. 
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example of high accretion rate objects (jBorosorul2002l ). but while they typically have steep 
X-ray slopes , their X-ray flux can range from X-ray bright to weak (e.g., = 0.7-2.2 
Leighlylll999l ). The T-Lboi/l-'Edd relation for the SDSS/XMM-Newton Quasar Survey will be 
discussed further in (Risaliti et al. 2009, in prep.). 

Inverse Compton scattering of UV photons from an accretion disk in a hot corona 
could explain why the X-ray spectrum pivots at low X-ray energies. In the corona model, 
the X-ray spectr um will change shape if the temperature (Tg) and/or optical depth (r) of 
the corona vary (IRybicki fc LightmanI Il979l ). since the output spectrum flattens as the y 
parameter increases. The y parameter is the average fractional energy gained by a photon; 
for a thermal, non-relativistic electron distribution. 



y 



max{T, r^) 



(2) 



For example, if the disk emission brightens, increasing the soft photon supply, the corona 
will increase radiative cooling to maintain temperature balance, producing a steeper X-ray 
spectrum. For a T ~ 10^-10^ K corona, opacity is not necessarily dependent on luminosity, 
so opacity variations c an result in a pivot in the 2-10 keV band without a large accompanying 
change in luminosity (IHaardt et al.l 119971 ). 



Because of the relation between the physical parameters and the output spectrum, it 
is possible to calculate TgT, though the d egenera cy is not b reakab le without knowing the 



cut-off of the high-energy spectrum. From iKrolik I (p. 227-8, Il999l ) 



-1.6 



where a is a coefficient dependent on geometry (a = 0.06 for slabs) and Ih/h is the compact- 
ness ratio, equivalent to the heating rate over the soft photon seed supply. Therefore, as F 
changes from 2.3 at low luminosities to 1.8 at high luminosities, TgX approximately doubles 
from 4.16 x 10® to 7.13 x 10® K. 

Several studies have found evidence of pivot points in individual objects. The previous 
paragraph describes the model of Mrk766, a NLSl, which was found to pivot near ~10 keV 



(IHaardt et al.lll997l ). In Cygnus X-1, a black hole X-ray binary, IZdziarski et al.l (120021 ) find a 
negative correlation between X-ray flux and F for the 0.15-12 keV band; no correlation for the 
20-100 keV band; an d a positive cor r elatio n for flux greater than 100 keV, implying a pivot 
point near ~50 keV. IZdziarski et al.l (120021 ) model this spectrum as a variable supply of soft 
seed photons irradiating a thermal plasma, where optical depth is constant. An increase in 
seed flux results in a decrease in the corona temperature to satisfy energy balance, resulting in 
a steeper spectrum with a pivot somewhere below 100 keV. Two mo del- independent analyses 
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of the spectral variability in Seyfert 1 galaxies also show evidence of pivot points. NGC 4051 
was fo und to have a pivot near 50 keV using correlations between the 2-5 keV and 7-15 keV 



f luxes ([Tavlor et al.l 20031 ). and NGC 5548, studied with an 8-day BeppoSAX observation 
( iNicastro et al.ll2000l ). can be entirely explained by a pivot in the medium energy band (~6 



keV). While these examples show that pivot points may exist in X-ray spectra, it is clear 
that further exploration must be done to provide a clear picture of the underlying physics. 



4.1.2. X-ray slope vs. optical color and aox (for RQ quasars only) 

X-ray slope (F) and optical color [A{g — i)] are correlated at 3.4(T, such that quasars 
with redder colors are more likely to have fiat X-ray slopes (Fig. 9a). This is likely due to 
undetected absorption, which is discussed further below. We exclude sources with redshifts 
greater than 2.3 from the F - A{g — i) correlation as Lya absorption artificially reddens 
high-redshift quasars. The linear bisector is: 



F = (2.04 ± 0.02) + (-0.45 ± 0.10)A(^ - i) 



(3) 



The T - A{g — i) correlation is consistent with previous fin dings of a correlation between F 
and the optical/UV s pectral slope (a,„,) (IKelly et al.l 120071 ) because color correlates tightly 
with the optical slope ( Richards et al. 20031 ) {A.{g — i) oc 0.5 a „„). The slope of th e correlation 
in this paper is in the same direction as, but steeper than in iKelly et al.l (120071 ). who find a 
slope = -0.25 ±0.07. 

X-ray slope (F) anticorrelates with aox at marginal (2.4cr) significance, such that X-ray 
faint quasars are more likely to have fiat X-ray slopes (Fig. 96). The WLS regression line is: 



F = (2.8 ±0.2) ± (0.6±0.1)ac 



(4) 



Th is anti-correl a tion b etween F and aox contrasts with the slightly positive correlation found 
bv I Green et al. I J2008h (overplotted as the almost horizontal dashed line in Fig. 86). The 
Green et al.l (120081 ) correlation may be affected by the inclusion of RL quasars, which lie at 
the X-ray bright end and typically have flatter X-ray slopes than RQ quasars. 

Completing the triangle of relations, a correlation is found between aox and A{g — i) 
at 4.7cr significance (Fig. 9c). Again, sources with redshifts greater than 2.3 are excluded 
in order to avoid contamination by the Lya forest. Since the aox-A{g — i) relation includes 
censored data, we use the EM method to give the best-fit line: 



On 



;-1.613 ± 0.009) ± (-0.16 ± 0.04)A(^ 



(5) 
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The EM algorithm assumes that all of the error lies in the dependent variable, so the algo- 
rithm reduces to an OLS Y vs. X fit for uncensored data. 

In an attempt to disentangle the relationships between F, A{g — i), and aox, we perform 
the Kendall partial correlation test for sources with S/N > 6 and z < 2.3, but each correlation 
is significant at the ~3cr level even when the third variable is taken into account. The 
dispersions in all three relations are large, so the sources are binned along the x-axis to 
show the relationships more clearly. The binning shows that the relations are dominated by 
outliers: optically red and X-ray weak/flat sources. 

Detected absorption cannot explain the observed correlations with F, since sources are 
only included if they do not prefer a model with intrinsic X-ray absorption. F is not correlated 
with Nh (Fig. 10a) for sources that prefer the APL model, indicating that once the X-ray 
slope is corrected for absorption, no intrinsic correlation between X-ray slope and absorption 
remains. For the same sources, aox and (Fig. 106) correlate at the 3.3cr level, which is 
surprising since the 2 keV flux should also be corrected for absorption. Since the detections 
trace the upper-limits, this correlation is likely not intrinsic, but is due to the survey flux 
limit. The optical color also depends on absorption (Fig. 10c) at the 2.6a level, as is expected, 
since X-ray absorbed sources are more likely to have redder colors due to dust-reddening. 

Undetected absorption is a possible explanation for all three relations, causing red opti- 
cal colors and X-ray weakness while flattening the X-ray spectral slope in low S/N spectra. 
A simple calculation shows that the amounts of absorption obtained from equations and 
(jl]) are consistent with observed gas-to-dust ratios observed for quasars. First, we assume 
that both relations are driven entirely by the effects of absorption (i.e. ignoring any possible 
effects due to intrinsic properties, such as black hole mass and accretion rate). We also 
assume that an unabsorbed quasar has a typical X-ray slope, (F = 1.9), zero absorption 
(Nh = and Eb-v = 0), and typical blue optical colors {A{g — i) = 0). If the X-ray slope 
changes by a given amount, we can calculate the cha nge in color by equation ([3]), which 



corresponds to a dust-reddening ([Richards et al.ll2003l ). which in turn leads to a change in 



optical luminosity. Since the X-ray slope also gives aox by equation (jl]), we can calculate the 
change in X-ray luminosity at 2 keV, which in turn gives the intrinsic X-ray absorption. So 
for an absorbed quasar with F = 1.4, the dust-reddening is E(B-V) = 0.096 from equation 
([3]), the gas column is ~ 6 x 10^^ cm~^ from equation (jlj), giving a gas-to-dust ratio ~ 
100 times the Galactic va lue. Quasars typically have gas-to-dust ratios in the range of 10 -100 



times the Galactic value (IMaccacaro et al.lll98ll : iMaiolino et al.ll200ll : IWilkes et al.ll2002i ). so 



relations ([3]) and (jl]) are consistent with intrinsic absorption that remains undetected in the 
X-ray spectra. 



Intrinsic absorption is less likely to remain undetected in spectra with high S/N. As a 
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second test, we restrict the correlation tests to those sources with X-ray S/N > 20. As a 
result, all three correlations disappear, which again supports undetected absorption as the 
cause of correlations that include low S/N sources. 



4.1.3. X-ray slope vs. redshift and optical luminosity 



Neither redshift nor 2500 A luminosity correlate si gnificantly with T (Ps > 10%, Figures 



11), confirming previ o us studies with smaller samples ( Page et al.l l2004l : iRisaliti et al 



2005 



Shemmer et al.l l2005l : IVignali et al.l l2005l : iKelly et al.l 120071 ) . Note however iBechtold et al. 



fl2003h 



who find that X-ray slopes are flatter at lower redshifts for a sample of 17 radio- 
quiet, high-redshift (3.7 < z < 6.3) quasars observed with Chandra. The SDSS/XMM sample 
(473 sources with X-ray spectra) is at least an order of magnitude larger than these previous 
samples, though it does not have homogenous spectral coverage for redshifts above 2.5. 



5. Conclusions and Future Work 

We have cross-correlated the DR5 SDSS Quasar Catalog with the XMM-Newton archive, 
creating a sample of 792 quasars with a detection rate of 87%. Almost 500 quasars have X-ray 
spectra, the largest sample available for analysis of optical/X-ray spectral correlations. We 
find that the X-ray photon index F correlates significantly with Ix, where X = 2, 4, 7, 10 and 
20 keV. Optical color and aox also correlate with F, but these correlations are hkely due to 
the effect undetected intrinsic absorption rather than intrinsic physical changes. The X-ray 
slope does not correlate significantly with redshift or optical luminosity. With a sample size 
at least an order of magnitude larger than previous studies, we confirm a highly significant 
correlation between aox and the monochromatic luminosity at 2500 A, and we also confirm 
that aox does not correlate significantly with redshift or optical color. 

Future studies of the sample will include: aox, F vs optical properties, and studies of 
sub-populations such as BALs and Type 2 quasars. Variability studies will also be pursued 
for 265 objects with multiple XMM observations. More complex spectral analysis on high 
X-ray S/N sources will include a thermal component for the soft excess, warm absorbers, 
and emission line detection. 
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6. Erratum: The Fifth Data Release Sloan Digital Sky 
Survey/XMM-NEWTON Quasar Survey (2009, ApJS, 183, 17); published 

Nov. 3, 2009 

We have discovered an error in column 9 of Table 2 in the original paper. This column 
reports the fit fiags, indicating which model a source prefers: a simple power-law (SPL, fiag 
= 1), a fixed power-law plus intrinsic absorption (FPL, fiag = 2), or an absorbed power- law 
(APL, fiag = 3). The table of the original paper mistakenly reports all fiag=3 sources as 
having fiag=2, and all fiag=2 sources as having fiag=l, so that only 32 sources prefer an 
absorption model. 



This preprint was prepared with the AAS lAT^TjX macros v5.2. 
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We have updated Table 2 to print out the correct fit fiags for each source, resulting in 
55 sources that prefer an absorption model. Since the fit fiags determine which numbers 
are reported for the remaining columns of Table 2, these numbers are updated as well. The 
abstract and text of the original paper report the correct number of absorbed sources, so the 
conclusions are unaffected. 

In addition, a minor rounding error was found in the SDSS names of some objects, and 
so we replace both Table 1 and 2 with corrected versions. Corrected machine-readable tables 
are available online from ApJS and will soon be available on the High Energy Astrophysics 
Science Archive Resource Center (HEASARC, http:/ /heasarc. gsfc. nasa.gov , ). 
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Table 1. SDSS quasars with archival XMM-Newton observations 





SDSS name 


XMM ID 


z 




(S/N)x^ 


Exp. 




Net 


Background 


RL 


BAL 














time'^ 




counts" 


counts^ 


fiags 


fiagi^ 


SDSS 


J125930.97+282705.5 


0204040101 


1.094 


0.92 


18.0 


221.0 


13.4 


469 


103. 








SDSS 


J130120. 13+282137.2 


0204040101 


1.369 


0.94 


45.8 


221.0 


13.2 


2780 


449. 


1 





SDSS 


J114856. 56+525425. 2 


0204260101 


1.633 


1.37 


23.0 


9.4 


5.8 


620 


52.1 


1 





SDSS 


J215419. 70-091753.6 


0204310101 


1.212 


3.71 


12.5 


81.3 


10.9 


218 


43.5 








SDSS 


J164221.22+390333.4 


0204340101 


1.713 


1.22 


11.3 


45.7 


11.9 


182 


39.1 








SDSS 


J021000.22-100354.2 


0204340201 


1.960 


2.20 


13.2 


31.6 


10.9 


243 


47.0 


1 





SDSS 


J021100.83-095138.4 


0204340201 


0.767 


2.17 


14.4 


31.6 


11.8 


268 


40.6 








SDSS 


J123508. 19+393020.0 


0204400101 


0.968 


1.49 


9.9 


65.5 


4.8 


153 


42.8 








SDSS 


J123527.36+392824.0 


0204400101 


2.158 


1.49 


20.1 


89.8 


6.0 


553 


104. 








SDSS 


J133812.97+391527.1 


0204651101 


0.439 


0.86 


18.7 


23.1 


8.0 


416 


38.4 









'^Galactic hydrogen column density (10^° cm^^) in the direction of the source. 
''X-ray signal-to-noise 
'^Exposure time (kiloseconds) 
'^Off-axis angle (arcminutes) 
"Net source counts 

^Background counts (counts in background region scaled to source extraction area) 

SRadio-loud flag is for RQ quasars, 1 for RL quasars, and 2 if the radio upper-limit is too high to determine whethe] 
source is RL. 

'^Broad-absorption line flag is for non-BAL quasars, 1 for BAL quasars 

Note. — Table [T] is presented in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown her( 
for guidance regarding its form and content. The upper limits are at the 1.6a level. 
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Table 2. X-ray Spectral Data of SDSS Quasars 



SDSSname 


to. 5-2 


12-10 


r b f b 
10.5-2 12-10 


Lo.5-2'^ L2_io'^ aox 


Fit 


po 








(obs.) 


(obs.) 


(rest) (rest) 


(rest) (rest) 











SDSS 


J125930.97+282705.5 


0.94 


5.20 


0.41 


2.46 


0.27 


1.61 


-1.77 


1 


n qq+0.13 


<3.10 


38.9/2, 


SDSS 


J130120.13+282137.2 


12.4 


20.5 


9.01 


17.0 


10.3 


19.3 


-1.32 


1 


1 70+0.07 


<0.05 


90.0/lS 


SDSS 


J114856.56+525425.2 


26.3 


59.9 


14.0 


39.5 


24.5 


69.4 


-1.60 


1 


1 58+°-i2 

l-3°-0.12 


<0.73 


20.3/21 


SDSS 


J215419. 70-091753.6 


1.51 


2.38 


0.93 


2.10 


0.79 


1.77 


-2.12 


1 


1 87+0.29 


<0.68 


5.0/12 


SDSS 


J164221. 22+390333.4 


3.09 


6.23 


1.76 


4.39 


3.48 


8.68 


-1.52 


1 


1.00_Q 24 


<0.44 


4.5/12 


SDSS 


J021000. 22-100354.2 


5.63 


12.7 


2.35 


8.04 


6.49 


22.1 


-1.37 


1 


1 60+°-i^ 

l-"^-0.19 


<1.38 


7.9/14 


SDSS 


J021100.83-095138.4 


5.72 


5.04 


5.72 


5.73 


1.56 


1.56 


-1.41 


1 


9 nq+0.24 
^•^"^-0.22 


<0.19 


5.6/ie 


SDSS 


J123508.19-f 393020.0 


1.68 


1.45 


1.74 


1.69 


0.84 


0.82 


-1.52 


1 


no+0.38 
■^■^'^-0.33 


<0.65 


5.1/lS 


SDSS 


J123527.36+392824.0 


2.57 


3.40 


1.57 


3.13 


5.49 


10.9 


-1.51 


1 


1 94+0.17 


<0.54 


20.7/2. 


SDSS 


J133812.97-f 391527.1 


9.13 


38.0 


6.69 


28.2 


0.47 


1.98 


-1.50 


1 


1-1 ' -0.16 


<0.08 


14.1/2: 



^Observed-frame X-ray flux in the observed band is given in units of 10^^^ ergs cm^^ s^^. 

'^Rest-frame X-ray flux in the soft and hard bands is given in units of 10^'^'' ergs cm"^ s~i. 

'^Rest-frame X-ray luminosity in the soft and hard bands is given in units of 10''^'* ergs s~i. 

flag indicating the X-ray fit. An undetected source is flagged as -1 and upper-limit flux values are listed. A deteci 
source with S/N < 6 is flagged as and detected flux values are listed. For sources with S/N > 6, three models can be appli 
If a single power-law model (SPL) is preferred, the flag = 1 and the SPL power-law parameters are listed, as well as the 9i 
upper-limit on intrinsic absorption from the intrinsically absorbed power-law (APL) model. If a flxed power-law (FPL) mo 
with variable Njy is preferred, the flag ~ 2. The best-fit slope from the SPL is listed, as well as the best-fit from the F 
model. If the APL model is preferred, the flag = 3, and the APL power-law and absorption parameters are listed. 

°Photon-index for the best-fit model. (If the APL model is preferred, then the photon index from that model is quot 
otherwise, the photon index is from the single power-law model.) 

^Intrinsic absorption or 90% upper-limit in units of 10^^ cm~^. 

^The value and degrees of freedom for the the best-fit model. 

Note. — Table [2] is presented in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown h 
for guidance regarding its form and content. Luminosities are computed using Ho = 70 km s-i Mpc-\ 17m = 0.3, and 
0.7. The upper limits are at the 1.6ct level. 
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Table 3. Weighted Averages of X-ray Spectral Quantities 







OLox 


N ^ 

-L> spec 


r 








All 


685 


-1.60 


473 


1.90±0.02 


0.40 


0.37 


0.15±0.03 


RQ+non-BAL 


589 


-1.61 


411 


1.91±0.08 


0.41 


0.38 


0.14±0.04 


RL+non-BAL 


62 


-1.46 


47 


1.85±0.03 


0.30 


0.26 


0.15±0.03 


BAL 


34 


-1.78 


15 


1.96±0.05 


0.37 


0.34 


2.3±0.6 



''Number of detected quasars in each sample (S/N > 2). 
•^Number of quasars in each sample with X-ray spectra (S/N > 6). 
■^Observed dispersion of X-ray slope. 
'^Intrinsic dispersion of X-ray slope 
^Intrinsic absorption in units of 10^^ cm~^. 



Table 4. F - Lx Correlations 



E (keV)'' 




Z-level^ 


slope 


intercept 


dispersion 


0.7 


0.0062 


2.73 


0.095±0.058 


-0.500±1.564 


0.36 


1.0 


0.1683 


0.16 


-0.015±0.054 


2.462±1.449 


0.37 


1.5 


0.0033 


1.77 


-0.106±0.051 


4.903±1.363 


0.37 


2.0 


0.0000 


2.93 


-0.162±0.048 


6.386±1.286 


0.37 


4.0 


0.0000 


5.61 


-0.265±0.041 


9.042±1.082 


0.35 


7.0 


0.0000 


7.50 


-0.316±0.035 


10.28±0.917 


0.33 


10.0 


0.0000 


8.60 


-0.335±0.032 


10.73±0.818 


0.31 


20.0 


0.0000 


9.61 


-0.365±0.027 


11.38±0.690 


0.29 



'^Energy at which monochromatic luminosity is taken. 

''Kendall's probability that correlation is due to chance. 

'^Significance level of Kendall rank correlation coefficient in units of Icr. 
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Fig. 1. — (a) Net counts vs. source extraction radius (arcseconds). (b) F vs. source 
extraction radius, (c) X-ray photon index (F) vs. off-axis angle (arcmin). (d) F vs. S/N. 
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Fig. 2. — Survey characteristics: (a) Exposure time histogram, where exposure times are 
summed over the up to three EPIC cameras. Exposure times do not include high-background 
intervals filtered out during data reduction, (b) S/N histogram, where 16 sources have X-ray 
S/N > 100. (c) Fraction of detected sources vs. redshift, and (d) redshift histogram. For (c) 
and (d), the open, dotted-hne histogram is for all SDSS-selected quasars, the open, solid-hne 
histogram is for all detected sources, the hatched histogram is for sources with sources with 
X-ray S/N > 6, and the double-hatched histogram is for sources with X-ray S/N > 10. 
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Table 5. T - Lx Correlations with redshift cut-offs 



E (keV)^ 


P b 




slope 


intercept 


dispersion 


0.7 


0.0112 


2 Ma 


0.583±0.285 


-13.29±7.587 


0.46 


1.0 


0.0386 


1.43a 


0.196±0.129 


3.077±3.447 


0.44 


1.5 


0.0024 


2.07(7 


-0.120±0.058 


5.270±1.363 


0.38 


2.0 


0.0000 


3.03cr 


-0.162±0.053 


6.800±1.417 


0.38 


4.0 


0.0000 


5.64(7 


-0.178±0.044 


9.596±1.176 


0.35 


7.0 


0.0000 


7.49cr 


-0.337±0.038 


10.84±0.986 


0.33 


10.0 


0.0000 


8.60(7 


-0.355±0.034 


11.25±0.874 


0.31 


20.0 


0.0000 


7.03(7 


-0.293±0.129 


9.560±3.342 


0.22 



'^Energy at which monochromatic luminosity is taken. 
''Kendall's probability that correlation is due to chance. 
■^Significance level of Kendall rank correlation coefiicient. 



100 



80 - 




xl 

Fig. 3. — xl histogram. The open histogram represents sources that prefer a power-law 
model. The solid histogram represents sources that prefer an absorbed power-law model. 
Five sources have > 2 and are not included in the plot. 
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^ 

-15 -14 -1,5 -12 

log F2-iOkev ergs cm"" s"' 



Fig. 4. — X-ray and optical sensitivity: (a) Optical 2500 A luminosity vs. redshift and (b) 
X-ray 2-10 keV luminosity vs. redshift. Flux limits are plotte d in both plots. For (a), i 
band magnitude hmits (dotted lines) were set during selection (IRichards et al.ll2002l ). For 
(b), t he flux limit (log F2- ^nkev = -14.0 ergs cm~^ s~^, dotted line) is for an observation of 
20 ks (IWatson et al.l 120011 ). (c) The observed- frame F2_iofcey distribution. 
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Fig. 5. — (a) hkev vs. /25oo- The dotted line shows hkev c<: ^2500, constrained to pass through 
the average values of hkev and /25oo- (b) The aox — ^2500 relation for the SDSS/XMM quasar 
sample. The linear fit to the SDSS/XMM sample is plotted as a solid line, and the Steffen 
et al. line is overp lotted as a dotted hne for comparison, (c) vs. redshift. Red circles 
represent detected sources with X-ray S/N < 6, and blue circles represent sources with X-ray 
S/N > 6. BALs and RL quasars were removed from all three plots. 



-2.5 -2.0 -1.5 -1.0 



1 



2 

r 



3 



Fig. 6. — (a) Q^ox distribution for all sources, (b) F distribution for sources with S/N > 6. 
Open histograms represent RQ, non-BAL quasars and have been scaled down by factors of 3 
and 4, respectively. Histograms with vertical hatching represent RL quasars, and histograms 
with diagonal hatching represent BAL quasars. The space above the upper-limit arrows 
indicates the number of RQ+non-BAL upper-limits in a given bar. 
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Fig. 7. — TOP: X-ray photon index (F) vs. monochromatic X-ray luminosity a t rest-frame 
(a) 0.7 keV, (b) 2 keV, (c) 10 keV, and (d) 20 keV. In (b), the lCreen et"aD JpOOsh c orrelation 
is plotted for comparison as a dashed line, and the quasars from Dai et al] ~f 2004 ) are over- 
plotted as open black squares. Only RQ, non-BAL quasars are plotted (blue circles) The 
weighted least squares regression line is calculated using RQ, non-BAL quasars, and is plotted 
as a solid line, with errors shown as dotted lines. BOTTOM: The weighted mean F values 
are given for bins of width Alog L2-iofcey = 0.5. Error lines mark the la dispersions in both 
axes. Blue lines represent RQ+non-BAL quasars, while green lines represent RL-|-non-BAL 
quasars. 
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Fig. 8. — The slope of the T-Lx correlation vs. energy E with la errors on the slope for 
(a) all sources and (b) sources where the rest-frame energy lies in the observed range. The 
slopes of the T-Lx correlations a re obtained for E = 0.7, 1, 1.5, 2, 4, 7, 10, and 20 keV, and 
are marked as .solid nirrle.s. The Oreen et a,l. mrrelatinn .slnne is marked as an nnen 
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Fig. 9. — TOP: (a) T vs. aox, with the T — correlation found by I Green et ahl (120081 ) 
plotted for comparison as a dashed line, (b) F vs. A{g — i), the relative {g — i) color, for 
sources with z < 2.3. Symbols for (a) and (b) are as in Figure U\ though the correlations 
are determined only for RQ, non-BAL and non-absorbed quasars. Weighted least square 
regressions are plotted as solid lines, with 3a errors shown as dotted lines, (c) aox vs. 
relative {g — i) color for sources with z < 2.3. Symbols for (c) are as in Figure O We plot 
the OLS regression as a solid line, with la errors shown as dotted lines. BOTTOM: The 
weighted mean F values are given for bin widths as follows: (a) Auox = 0.2, (b) A[A{g — i)] 
= 0.3, and (c) A[A{g — i)] = 0.4. Error lines mark the la dispersions in both axes. Blue lines 
represent RQ+non-BAL quasars, while green lines represent RL+non-BAL quasars. Figure 
(c) does not include BAL or RL quasars. 
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Fig. 10. — (a) Intrinsic absorption (Nh) vs. T. (b) Nj^ vs. agx- (c) vs. A{g — i). 
The solid circles in (a) and (b) represent sources preferring the APL model, while in (c) the 
solid circles represent sources preferring either absorption model (FPL or APL). The arrows 
represent sources with only a 90% upper-limit on Nj^. Nj^ is given in log units of 10^^ cm~^. 
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Fig. 11. — X-ray photon index (F) vs. (a) optical 2500 A monochromatic luminosity and 
(b) redshift. Blue circles represent RQ, non-BAL quasars, green stars represent RL quasars, 
and red triangles represent BAL quasars. Filled symbols are given for those sources that 
prefer an absorbed power-law with an F-test probability P > 0.95. The bottom plots show 
the weighted mean F values for bins of width (a) Alog L2500 = 0.5 and (b) Az = 1. 



